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Abstract 

A new dynamical model is developed to describe the whole process of surrogate reactions; trans- 
fer of several nucleons at an initial stage, thermal equilibration of residues leading to washing out 
of shell effects and decay of populated compound nuclei are treated in a unified framework. Multi- 
dimensional Langevin equations are employed to describe time-evolution of collective coordinates 
with a time-dependent potential energy surface corresponding to different stages of surrogate reac- 
tions. The new model is capable of calculating spin distributions of the compound nuclei, one of the 
most important quantity in the surrogate technique. Furthermore, various observables of surrogate 
reactions can be calculated, e.g., energy and angular distribution of ejectile, and mass distributions 
of fission fragments. These features are important to assess validity of the proposed model itself, 
to understand mechanisms of the surrogate reactions and to determine unknown parameters of the 
model. It is found that spin distributions of compound nuclei produced in 18 0+ 238 U — )> 16 O+ 240 *U 
and 18 0+ 236 U — ^ 16 0+ 238 *U reactions are equivalent and much less than lO/i, therefore satisfy 
conditions proposed by Chiba and Iwamoto (PRC 81, 044604(2010)) if they are used as a pair in 
the surrogate ratio method. 

PACS numbers: 24.87.+y, 24.10.-i, 24.60.Ky, 27.90.+b 
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I. INTRODUCTION 



Neutron-induced cross section data of unstable nuclei are systematically required to design 
next-generation nuclear facilities such as high burn-up fast breeder reactors or accelerator- 
driven systems for transmutation of nuclear wastes [1], Such data are also important 
to understand origin of elements, namely, the s- and r-process nucleosynthesis (see, e.g., 
Refs. [2I, S]). However, it is not usually possible to measure these cross sections directly by 
using neutrons due to difficulty in preparing samples. It is well known that the compound 
reaction process is the dominant mechanism in the energy region of our interest. There- 
fore, various experimental methods to measure the "direct" neutron capture components 
are not be applicable to determine them. Instead, other methods are needed, and one of 
the promising method is the surrogate reaction approach [^-[25!]. In this method, (multi) 
nucleon transfer reactions with an experimentally accessible combination of projectile and 
target are employed to create the same compound nucleus as the desired neutron reaction, 
and decay branching ratios to specific channels, normally capture and/or fission, are de- 
termined. However, branching ratios are sensitive to the spin and parity of the compound 
state, while the spin-parity (J 71 ") distributions of populated nuclei are probably different for 
the neutron-induced and surrogate reactions. It is well known that if the spin is different 
even just 1 unit, the capture branching ratio is totally different for the energy region of our 
interest. Therefore, validity of the surrogate method depends on how the difference of the 
spin-parity distributions is comprehended and compensated properly. 

Recently, the surrogate ratio method (SRM) is discussed by Chiba and Iwamoto. It was 
found that SRM works to a certain accuracy if (1) there exist two surrogate reactions whose 
spin-parity distributions of decaying nuclei are equivalent, (2) difference of representative 
spin values between the neutron-induced and surrogate reactions is not much larger than 10 
h, under a condition that (3) weak Weisskopf-Ewing condition, namely, J^-by-J^ conver- 
gence of the branching ratio, is realized [26|]. They form a set of sufficient conditions for the 
SRM to work. It is important to notice that the J 71 " distribution may be even different for 
the neutron-induced and surrogate reactions if these conditions are fulfilled. 

Discovery of the above conditions is a great advancement for the whole surrogate tech- 
nique. Therefore we need further investigation to verify that the above conditions, especially 
(1) and (2) which were just assumed in Ref. [26|], are really satisfied in certain surrogate 
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reactions. It implies that mechanisms of the surrogate reactions to be understood well. For 
that aim, it is indispensable to establish a theoretical model to describe the whole process of 
the surrogate reactions, namely, nucleon transfer and decay of populated compound nucleus. 
In this work, we propose a first version of our model to describe the surrogate reactions 

n 

based on a theory proposed originally by Zagrebaev and Greiner [23] . This model, called 
a unified model, can treat the whole reaction processes in heavy- and superheavy-mass 
regions, which has been applied to several types of reactions [27H30J. The name of unified 
model implies an unified dynamical approach and unified multidimensional potential energy. 
Time-evolution of the system is described by a trajectory calculation on the time-dependent 
unified potential energy surface using the Langevin equation. Then, we treat a two neutron 
transfer reaction; 18 0+ 238 U — >> 16 O+ 240 U, which is planned to be employed at Japan Atomic 
Energy Agency (JAEA) as an example of application of the new model to the surrogate 
reaction. By using the new model, we can obtain various quantities which can be compared 
with experimental data directly and we can evaluate our theory and determine unknown 
parameters in the model. 

The purpose of this paper is to explain the new model, calculate a fission fragment 
mass distribution (FFMD) for a reaction 18 0+ 238 U to calibrate the model parameter, and 
calculate spin distributions of compound nuclei for the 18 0+ A U — ^ 16 0+ A+2 *U system, where 
A =236 and 238, to see if the conditions proposed by Chiba and Iwamoto are satisfied or not. 
In section 2, we explain our theoretical framework. The calculation results are presented in 
section 3. In section 4, we present a summary of this study and further discussion. 



II. DYNAMICAL MODEL 



Overview of the model 



The surrogate reactions consist of 2 stages; an initial nucleon transfer process and de- 
cay of populated compound nuclei, which have quite different nature to each other. The 
Hauser-Freshbach (HF) theory 0, 32] has been applied to describe the latter part of the 
surrogate reactions. By the HF theory, we are able to calculate decay branching ratios to 
specific channels (capture or fission) with arbitrary spin-parity distributions of compound 
states. However, we cannot predict the spin distribution produced by the initial stage of the 
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surrogate reaction nor the FFMD with the HF theory. For that reason, we need to describe 
the whole reaction process consistently, i.e., starting from the transfer of several nucleons 
and the decay of the compound nucleus leading to fission successively. Here, we employ a 
dynamical model, the unified model, to the surrogate reaction. 

The unified model was proposed by Zagrebaev and Greiner and was applied to several 



types of reactions 
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-30| induced by heavy ions. An unified dynamical approach and unified 
multidimensional potential energy are employed in this model, which are the origin of the 
name of this theory. To apply this model to the surrogate reaction, we extend this model 
and introduce new procedures. As explained above, the surrogate reactions consist of two 
processes; the transfer reaction process between a two-body system and decay of the pop- 
ulated compound nuclei (one-body system), for which the mass of total system is different 
very much to each other. Therefore, it is indispensable to connect such different systems 
to treat surrogate reactions consistently. Evolution of the mass-asymmetry parameter is 
described by multi-dimensional Langevin equations without (with) the inertia parameter 
before (after) the window of the colliding nuclei opens sufficiently. We modify the original 
unified model [22j also to take account of temperature dependence of the shell correction 
energy of the potential energy surface. 

Firstly, we treat the transfer reaction process within the framework of the unified model. 
Then, we treat the decay of the compound nuclei with an initial condition populated by the 
former reaction process. We perform a trajectory calculation on a time-dependent potential 
energy surface corresponding to different stages of the surrogate reaction. A dynamical 
calculation is carried out in terms of the multi-dimensional Langevin equation based on the 
fluctuation-dissipation theorem. By this procedure, we can describe trajectories moving on 
the potential energy surface including the nucleon transfer reaction. 

We consider in this paper a two neutron transfer reaction 18 0+ 238 U — >> 16 O+ 240 *U as an 
example. In the transfer reaction process, we use the potential energy of 256 Fm as the total 
system. After the production of the compound nucleus by the transfer reaction, we then 
treat the decay of the compound nucleus 240 U. In other words, the potential energy surface 
switches from that of 256 Fm to that of 240 U as the surrogate reactions proceed from the 
transfer process to the decay stage. 

A schematic picture of our model is presented in Fig. [H It shows the potential energy sur- 
faces used in the trajectory calculation from the transfer process to the decay of compound 



nuclei in the reaction 18 0+ 238 U — )> 16 O+ 240 U. The transfer reaction process is presented in 
the upper left panel of Fig. 1 which shows a diabatic potential energy surface of 256 Fm in 
the z-ol (5 = 0) coordinate space. Meanings of these parameters and terms are explained 
in the next subsection. The thin arrows correspond to the entrance and the exit channels 
of two-neutron transfer process. In the calculation, it starts from an infinite distance be- 
tween the projectile and target, where actually the distance of 30 fm is used. Then, the 
calculation stops when the trajectory reaches at the distance of 25 fm between the both 
fragments. We select the two nucleon transfer among the all events; we select events in 
which the mass asymmetry parameter a changes from 0.859 corresponding to 18 0+ 238 U to 
0.875 corresponding to 16 O+ 240 U. 

Next, the decay process of compound nucleus is presented in the lower right panel of 
Fig. 1, which is a potential energy surface of 240 U with S = 0.2 (/3 2 ~ 0.2). The white lines 
denote mean fission paths. In the decay process, we start the trajectory calculation with the 
initial condition obtained in the transfer reaction process. As the initial condition, we use 
the deformation, momentum, angular momentum and excitation energy of the compound 
nucleus. The quantities which we can obtain with this calculation are the angular and energy 
distributions of ejectile, mass and total kinetic energy distributions of fission fragments and 
neutron multiplicities, and so on. Such quantities can be compared with experimental data 
directly, which allows us to determine unknown parameter in the model. In this work, we 
discuss a mass distribution of fission fragments and determine a value of one of the most 
uncertain parameter, the sliding friction (see below). In this way, we can evaluate and 
improve our model step-by-step by comparing model predictions with experimental data. 
More details of our model will be explained below. 

B. Potential energy surface 

The initial stage of the surrogate reaction consists of 2 parts; 1) a fast diabatic part in 
which the reaction proceeds too fast for nucleons to reconfigure their single-particle states 
so the system goes through the ground-state configurations of the target and projectile, and 
2) the system relaxes to the ground-state of the total composite system, which changes the 
potential energy surface to an adiabatic one. Therefore, we take into account time evolution 
of the potential energies from the diabatic one Vdidb(q) to the adiabatic one V a di a b(q), here 
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q denotes a set of collective coordinates representing the nuclear deformation. The diabatic 
potential is calculated by a folding procedure with effective nucleon-nucleon interaction 

3, 



28l . l33| . which is shown in the upper left part of Fig. [TJ We can see a "potential wall" 
in the overlap region of the colliding system, which corresponds to a hard-core representing 
incompressibility of nuclear matter. On the other hand, the adiabatic potential energy of 



the system is calculated using an extended two-center shell model 



331 ] . We then connect the 



diabatic and adiabatic potentials with a time-dependent weighting function as follows; 

V = V diab (q)f(t) + V adiab (q)[l - f(t)}, 

f(t) = exp (--) . (1) 



T j 

Here, t is the time of interaction and f(t) is the weighting function with the relaxation time 
r. We use a relaxation time r = 10 -21 sec, which was suggested in references j^]. It is 
empirically known that calculated results do not depend noticeably on the relaxation time. 
As the coordinates to express nuclear deformation, we use the two-center parametriza- 



tion 



36| and employ three parameters as follows: z G (distance between centers of two 
potentials), 8 (deformation of fragments), and a (mass asymmetry of the colliding nuclei); 
a = (Ai — A 2 )/(Ai + Ao), where A\ and A 2 denote the mass numbers of the target and the 



projectile, respectively (37j. Later on, A\ and A 2 are used to denote mass numbers of two 
fission fragments. The parameter 8 is defined as 8 = 3(a — b)/ (2a + 6), where a and b denote 
the half length of the axes of ellipse in the z and p direction, respectively as expressed in 
Fig. 1 in reference [35]. We assume that each fragment has the same deformations as a first 
step. Furthermore, we use scaling to save computation time and employ a coordinate z de- 
fined as z = Zq/(RcnB), where Ron denotes the radius of the spherical compound nucleus. 
The parameter B is defined as B = (3 + 8)/ (3 — 28). 

In the two-center parametrization, the neck parameter is denoted by e and is known to 
be different in the entrance and exit channels j^]. Therefore, we employ e — 1 for the 
entrance channel and e = 0.35 for the exit channel to describe a realistic nuclear shape. We 
introduce a time-dependent potential energy surface in terms of e using a relaxation time 
for e of r e = 10 -20 sec f^ . as follows; 

V a diab = V adiab (q, e = l)f e (t) + V adiab (q, e = 0.35) [1 - f € (t)], 

fe(t) = e W f~Y (2) 
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C. Dynamical equations 



We then perform trajectory calcu' 



using the Langevin equation 



ations on the time-dependent unified potential energy 
33 



The nucleon transfer for slightly separated nuclei is important in surrogate reactions. 
Such intermediate nucleon exchange plays an important role in fusion process at incident 
energies near and below the Coulomb barrier as well. We treat the nucleon transfer using 



the procedure described in reference [221, l28[ ; 



da 
~dt 



A 



CN 



-D%\a) + —^D%\a)T a (t), 



CN 



(3) 



This is a Langevin equation neglecting inertia mass for the mass asymmetry parameter. 
It expresses change of the asymmetry parameter a due to drift (first term on the right 
hand side) and diffusion (second term on the r.h.s.) processes. It is obtained by a certain 
approximation starting from the Master equation for transition of different particle-hole 
states. Such Master equation giving discrete change of nucleon numbers is transformed to 
an equation for a continuous variable a via Fokker-Planck equation to the Langevin equation 
shown above [271. 1281]. 



After the window of the touching nuclei opens sufficiently (hereafter "the mono-nucleus 
state"), the treatment of the evolution of the mass- asymmetric parameter a switches from 
eq. (EJ) to the Langevin equations with the procedure described in reference js^j. The 



multidimensional Langevin equations 27, [jol, Q are now unified as 

d>Qi ( —\\ 

dt v J v Fv 

dpi dV 1 d , _^ , _ n 

+9ijRj(t), 

de_ _ ^ 

(it /X#i? 2 ' 
dt ' 

dt dV ( £ L x L 2 \ 

It = ~ ltan \J^R ~ ~ % a y R 

dL r dV ( i Li L 2 \ 



~ a l9tanRtan(t) , 

rfL 2 ay ( t l x l 2 



dt d(f2 \^rR ^1 ^2 

-Q>29tanRtan{t), (4) 

where a summation over repeated indices is assumed. The collective coordinates % stands 
for z, 5 and a. The symbol denotes momentum conjugate to and 1/ is the multi- 
dimensional potential energy. Definition of other parameters is given in Fig. 2: The symbols 
6 and £ are the relative orientation of nuclei and relative angular momentum, respectively, 
and (fi and denote the angles of rotation of the nuclei in the reaction plane (their moments 
of inertia and angular momenta are and L^ 2 , respectively), = R/2 ± (Ri — ^)/2 
are the distances from the centers of the fragments up to the middle point between nuclear 
surfaces, and Ri^ are the nuclear radii. The symbol R is the distance between the nuclear 
centers. The total angular momentum L = £ + Li + L 2 is conserved. The symbol jiR denotes 
the reduced mass, and j tan is the friction force in the tangential direction of colliding nuclei, 
here we call it as the sliding friction. 

The symbols and 7^ stand for elements of the shape-dependent collective inertia 
and friction tensors, respectively. For separated nuclei, we use the reduced mass and the 
phenomenological friction forces with the Woods-Saxon radial form factor as described in 



reference 



281 ]. We switch the phenomenological friction to the friction for mono- nuclear 



system using a smoothing function 



281 ]. For the mono- nuclear system, the wall-and- 



window one-body dissipation is adopted for the friction tensor, and a hydro dynamical 



ertia tensor is adopted in the Werner- Wheeler approximation for the velocity field [39l-l4ll] 



m- 



The normalized random force Ri(t) is assumed to be of white noise, i.e., (Ri(t))=0 and 
(Ri(ti)Rj(h)) = 2SijS(ti — £2)- Strength of the random force g^ is given by Einstein re- 
lation; jijT = ^2 k gijgjk, where T is the temperature of the compound nucleus calculated 
from the intrinsic energy of the composite system. 
The adiabatic potential energy is defined as 

V adiab (q, L, T) = V LD (q) + 1} + V SH (q, T), (5) 

V LD (q) = E s (q) + E c (q), (6) 

V SH (q,T) = E° shell (lMT), (7) 

$(T) = exp (-!p\ , (8) 



, Ed, 

where I(q) stands for the moment of inertia of a rigid body with deformation q : Vld and 
Vsh are the potential energy of the finite-range liquid drop model and the shell correction 
energy taking into account the temperature dependence, respectively. The symbol E® hell 
denotes the shell correction energy at T = 0. The temperature dependent factor $(T) is 
discussed in reference [421 ]. where E* denotes the excitation energy of the compound nucleus. 
The shell damping energy E^ is chosen as 20 MeV, which is given by Ignatyuk et al. [431 ]. 

The symbols Es and Ec denote a generalized surface energy [44] and Coulomb energy, 
respectively. The centrifugal energy arising from the angular momentum L of the rigid body 
is also considered. The intrinsic energy of the composite system E int is calculated for each 
trajectory as 

E^t = E*~\ (m" 1 ) . . p iPj - V(q, L, T). (9) 

Here, E* is given by E* = E cm — Q ) where Q and E cm denote the Q- value of the reaction 
and the incident energy in the center-of-mass frame, respectively. Each trajectory starts 
from a sufficiently large distance between both nuclei [301 ]. 



D. Computation 



Due to difference of the initial impact parameters (or the different initial relative angular 
momenta), various kind of reaction processes can occur. Moreover, even the trajectories start 
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with the same initial impact parameter, reactions proceed in quite different ways due to the 
random force in the dynamical equation, which finally leads to different reaction channels. 
By choosing various impact parameters randomly and give proper weights, whole processes 
of reactions are described by the present model; the elastic and inelastic scattering, deep 
inelastic collision, quasi-fission, fusion-fission process, and a few nucleon transfer process (the 
surrogate reaction). They are treated simultaneously by the model. This is a big advantage 
of the present approach, since these reactions correlate, thus they can give information to 
each other. An example is a determination of the unknown parameter 7 tan through FFMD 
of fusion-fission like process as will be discussed below. 



III. RESULTS AND DISCUSSION 



The unified model, on which the present model is based, has been applied to several types 



27- 



30| . In our previous study, 



of reactions and succeeded to describe the experimental data 
we precisely investigated the incident energy dependence of mass distribution of fission 
fragments in the reactions 36 S+ 238 U and 30 Si+ 238 U |30|, 



45, 



46( | . The calculation results 



reproduced the experimental data well and clarified the origin of the fine structure of the 
mass distribution of fission fragments at the low incident energy. 

Here, we apply the present model to the surrogate reaction. To evaluate and clarify 
the model, we focus on the mass distribution of the fission fragments and compare calcu- 
lated results with experimental data. Furthermore, we investigate spin distributions of the 
compound nucleus populated by transfer reactions and discuss the validity condition of the 



SRM 



26|]. We choose a system of 18 0+ A U -V 6 0+ A+2 *U reaction where A =236 or/and 238. 



In the Langevin calculation, the sliding friction is mainly responsible for the dissipation 
of the angular momentum |48|, |49|, though its value is uncertain. In the present work, 
we treat the sliding friction as a parameter of the model and investigate dependence of 
the calculation upon this parameter. Measured FFMD data in the reaction 18 0+ 238 U at 
E c .m = 133.5 MeV is shown by dots in Fig. [31 It includes all the fission fragments occurring 
in the above reaction. The experimental set-up and the data analysis are nearly the same as 
the references (45, 47]. In the experiment, both fission fragments were detected in coincidence 
by using position-sensitive multiwire proportional counters (MWPCs). The difference of the 



setup from the references 



471 ] was the angles of the detector positions that MWPC1 and 
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MWPC2 were located at —90° and +90° with respect to the beam direction. 

In Fig. [3l calculated data with 7 tan = 1, 5, 10 and 20 xlO -22 MeV s fm~ 2 are denoted by 
the black, red, green and dark blue histograms, respectively. It is clearly seen that the result 
with ^ tan = 5 reproduces the experimental data very well. With larger sliding friction, the 
variance of results becomes smaller. We can fix values of unknown parameters in this way. 

Then, calculations for the 2 neutron transfer reaction are carried out without any ad- 
justable parameters. From such calculations, we can determine spin distributions of the 
surrogate reactions 18 0+ 236 U i6Q+ 238 *U and 18 0+ 238 U -> 16 O+ 240 *U, and verify if the 
2 assumptions proposed by Chiba and Iwamotoj^, namely, (1) there exist two surrogate 
reactions whose spin-parity distributions of decaying nuclei are almost equivalent, and (2) 
difference of representative spin values between the neutron-induced and surrogate reactions 
is not much larger than 10 H : are really satisfied or not. 

Figure H] shows calculated spin distributions of compound nucleus 240 U populated in the 
reaction 18 0+ 238 U — )> 16 O+ 240 U at an incident energy of E c , m = 160 MeV, which is planned 
to be carried out at JAEA. Results with various values of the sliding friction are shown. We 
can see that majority of the spin of compound nucleus is much less than 10 h for each value 
of the sliding friction although they diverge depending on ^ tan . Therefore, it is important for 
the model to have a capability to determine values of unknown parameters as shown above. 
Our model is particularly powerful since this parameter is determined by using observables 
corresponding to other reaction channels, which can be treated simultaneously with the 
surrogate reactions of our interest. Figure [5] shows spin distributions of compound nuclei 
240 U and 238 U in the transfer reactions 18 0+ 238 U -^ 16 O+ 240 U and 18 0+ 236 U -V 6 0+ 238 U, 
respectively with the sliding friction ^ tan = 5 x 10 -22 [MeV s fm -2 ]. These distributions 
should be interpreted as a semi-classical estimate of a spin distribution corresponding to 
excitation of rotational motion due to angular momentum transfer occurring in these reac- 
tions. It is easily noticed that the spin distributions of the compound nuclei populated by 
the two reactions are almost equivalent. These results suggest that the assumptions (1) and 
(2) shown above for the SRM to work[26] are proved to be correct within this model. In 
conjunction with the weak Weisskopf-Ewing condition proposed in Ref. [26] (which can be 
verified by Hauser-Feshbach theory), the present result suggests that 18 0+ 238 U — >> 16 O+ 240 U 
and 18 0+ 236 U ^ 16 0+ 238 U reactions can be employed as a pair in the SRM. 

We propose that this model as a powerful and useful tool to describe the surrogate 
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reaction process, even though it is a semi-classical model. To describe the transfer reaction 
process more accurately, we may have to consider quantum effects precisely. Quantum 
mechanical models such as DWBA or CDCC [5C| would be more appropriate to describe 
the nucleon-transfer part of the surrogate reactions. It will be absolutely necessary if we 
use light-ion projectiles. We can use these sophisticated models if they are available and 
connect the populated spin-distribution to the later part of the present model. However, 
it is difficult to treat transition probabilities to continuous levels quantum mechanically as 
realized in surrogate reactions. As the first step, therefore, we try to understand the gross 
feature of the surrogate reaction and analyze the reaction mechanism using the present 
model in this paper. Especially, it is known that the dynamical model is useful to discuss 



the mass distribution of fission fragments 



51 



521 ] , an important observable of the surrogate 



reactions which contains information on the populated compound nuclei such as the spin 
distribution If 



IV. SUMMARY 

We propose a first version of a unified dynamical theory to describe the whole process 
of surrogate reactions; the nucleon transfer, thermalization and the decay of the popu- 
lated compound nuclei. To realize it, we introduced new procedures to the unified theory 
of Zagrebaev and Greiner, namely, switching of the potential energy surfaces having very 
different mass numbers, Langevin equations depending on different stages of the reaction 
and a temperature-dependent shell correction energy. Trajectory calculations in terms of 
the Langevin equations are employed on a time-dependent potential energy surface corre- 
sponding to different stages of the surrogate reactions. After the transfer process, decay of 
the populated compound nucleus is calculated with the initial condition obtained from the 
preceding transfer process. This model can yield many observables which can be compared 
with experimental data directly. 

As an example of the application of the present model to surrogate reactions, we con- 
sidered a two nucleon transfer reaction; 18 0+ 238 U — )> 16 O+ 240 U, which is planned to be per- 
formed at JAEA. We treated the sliding friction as a parameter of the model and discussed 
the dependence of the calculation results upon the sliding friction. Then, we discussed the 
validity condition of the surrogate ratio method (SRM). We calculated the spin distribution 
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of the compound nuclei with several sliding frictions for the compound nucleus 240 U in the 
reaction 18 0+ 238 U -^ 16 O+ 240 U at the incident energy of E c , m , = 160 MeV. The calculation 
results showed that the spin of compound nucleus was less than 10h for each value of the 
sliding friction. Finally, we discussed spin distributions of compound nucleus 240 U and 238 U 
in the transfer reactions 18 0+ 238 U -V 6 O+ 240 U and 18 0+ 236 U -V 6 0+ 238 U, respectively. It 
was found that the spin distributions of decaying nuclei populated by the two reactions are 
almost equivalent. Therefore it is concluded that if these reactions are used as a pair in 
the SRM, they would yield the correct neutron cross sect ions [26]. These calculation results 
suggested validity of the SRM within this model. 

In the present model, however, nuclei were treated as nuclear matter, and a semi-classical 
approach was employed except for the fact that we took into account the shell correction 
energy on the potential energy surface. Such semi-classical model may be too simple, and 
we may have to consider quantum effects in order to describe the reactions more accurately. 
Nevertheless, the present model is flexible enough to take account of results of more elab- 
orated models. We therefore consider that the present model is capable enough, as a first 
step, to understood gross features of the surrogate reactions which itself is already quite 
complicated. 

As further studies, we improve the model in the nucleon transfer part by taking into 
account the quantum effect more precisely. After experiments of the surrogate reaction 
at JAEA are performed, we can compare model predictions with experimental data, e.g., 
distributions of emission angle and energy loss of ejectile, mass, charge and total kinetic 
energy distributions of fission fragments from various exclusive fission channels, and the 
model will be upgraded successively. 
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FIG. 1. (Color online) Schematic picture of the calculation. For the surrogate reaction 18 0+ 238 U 
— ^ 16 O+ 240 U, the potential energy surfaces from the transfer reaction to the decay process of com- 
pound nuclei are presented. The transfer reaction is shown in the left panel which is the diabatic 
potential energy surface of 256 Fm in the z — a(S = 0) coordinate space. The decay process of 
compound nucleus is presented in the right panel, which is the adiabatic potential energy surface 
of 240 U with 6 = 0.2 (/3 2 - 0.2) . 
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FIG. 2. (Color online) Definition of parameters used in the model. 
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FIG. 3. (Color online) Fragment mass distribution obtained in the reaction 18 0+ 238 U at an incident 
energy of E CmTrim = 133.5 MeV. Experimental data and calculation results are denoted by circles and 
histograms, respectively. Calculations are shown with sliding frictions ^y t an — 1, 5, 10 and 20 x 10 -22 
[MeV s fm -2 ], which are multiplied by the factor such that the total cross section agree with the 
experimental value to compare the shape of the mass distribution with the experiment. 
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FIG. 4. (Color online) Spin distribution of compound nucleus 240 U in the reaction 18 0+ 238 U 
^16q + 240 U at the i nc ident energy of E CmTrim — 160 MeV for several sliding frictions. The black, 
red, blue and green lines denote for ^y tan = 1, 5, 10 and 20 xlO -22 [MeV s fm -2 ], respectively.. 
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FIG. 5. (Color online) Spin distribution of compound nuclei 240 U and 238 U in the transfer reactions 
18q + 238tj _^16 O+ 240tj and 18q + 236tj _ )> 16q + 238 U ^ at the incic ient energy of E c . m . = 160 MeV, 

respectively. Sliding friction 7 tan = 5x 10 -22 [MeV s fm -2 ] is used. 
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